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Currently the theoretical uncertainty limits the interpretation of the atomic parity non- 
conservation (PNC) measurements. We calculate the PNC 5s - 6s electric dipole transition am- 
plitude in rubidium and demonstrate that rubidium is a good candidate to search for new physics 
beyond the standard model since accuracy of the atomic calculations in rubidium can be higher 
than in cesium. PNC in cesium is currently the best low-energy test of the standard model, there- 
fore, similar measurements for rubidium present a good option for further progress in the field. We 
also calculate nuclear spin-dependent part of the parity non-conserving (PNC) amplitude which is 
needed for the extraction of the nuclear anapole moment from the PNC measurements. 

PACS numbers: 11.30.Er, 31.15.A- 
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INTRODUCTION 



The study of parity non-conservation in atoms is a low- 
energy, relatively inexpensive alternative to high-energy 
searches for new physics beyond the standard model (see, 
e.g. [U-IH). For example, parity non-conservation in ce- 
sium is currently the best low-energy test of the elec- 
troweak theory Q], Q . It is due to high accuracy of the 
measurements [j] and its interpretation [B| (see also @- 
8]). The uncertainty of the measurements is 0.35% 
while the uncertainty of the calculations is on the level 
of 0.4% - 0.5% MB- This means that the interpreta- 
tion of the measurements is limited by the accuracy of 
atomic calculations. The situation is similar for other 
atoms. For example, accuracy of the PNC measurements 
for thallium is 1% [9j while accuracy of the calculations is 
2.5-3% [Id, EH. It is believed that a good option for fur- 
ther progress may come with the PNC measurements for 
atoms or ions with electron structure similar to cesium 
but with higher nuclear charge Z. Higher Z would lead 
to larger PNC effect and would probably lead to better 
accuracy in the measurements. The PNC measurements 
have been considered for the Ba + ion 12, Il3j and are 



under progress for the Fr atom [l4| and Ra + ion [If 
However, the accuracy of the calculations for these sys- 
tems is unlikely to be better than for cesium. Just on 
the contrary, higher Z means larger relativistic effects 
such as Breit and quantum electrodynamics (QED) cor- 
rections, larger uncertainty due to the neutron skin effect, 
etc. This would most likely lead to poorer accuracy of 
the calculations. Since the accuracy of the calculations 
is a limiting factor even for cesium, it does make sense 
in our view to look in opposite direction and to consider 
PNC in rubidium. Rubidium is an alkali atom next to 
cesium but with smaller Z. The same calculations as for 
cesium would lead to better accuracy of the results for 
rubidium while the PNC amplitude in rubidium is only 
seven times smaller than in cesium. Depending on the 
accuracy of the measurements which can be achieved for 



rubidium, the study of the PNC in this atom might be a 
good alternative for further progress in the area. 

Rubidium was considered for anapole moment mea- 
surements in Ref. L 16j]. Corresponding atomic calcula- 



tions were reported in Ref. [13, LL8J . The calculations of 
the spin-independent PNC amplitude of the 5s - 6s elec- 
tric dipole transition in Rb were performed in our early 
work 1 10] with 2% accuracy. Only correlation corrections 
were considered while Breit, QED and other small cor- 
rections were ignored. In this paper we perform a more 
detailed analysis of the PNC amplitude in Rb. This in- 
cludes more accurate treatment of the correlations and 
a detailed consideration of the Breit, QED and neutron 
skin corrections. We demonstrate that the analysis of the 
PNC effect in rubidium can be more reliable and accurate 
than in cesium. This should encourage experimentalists 
to consider the PNC measurements in rubidium. 



CALCULATIONS 

The Hamiltonian describing parity-nonconserving 
electron-nucleus interaction can be written as a sum of 
the nuclear-spin-independent (SI) and the nuclear-spin- 
dependent (SD) parts (we use atomic units: H = |e| = 
m e = 1): 
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where Gf ~ 2.2225 x 10~ 14 a.u. is the Fermi constant 
of the weak interaction, Qw is the nuclear weak charge, 

^ *T ] and 75 are the Dirac matrices, I is the 

<T J 

nuclear spin, and p(r) is the nuclear density normalized 
to 1. 

Within the standard model the weak nuclear charge 
Qw is given by 



-0.98777V + 0.0716Z. 



(2) 
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TABLE I: Ionization energies (in cm 1 ) and hyperfine struc- 
ture constants A (in MHz) for low states of 85 Rb. 



State 


Energies 


[cm- 1 ] 


A [MHz] 






Exper." 


Calc. 


Exper. 


Calc. 


0*1/2 


ooOy 1 


OoDDD 


lun.y 


1U1D 


5^1/2 


21113 


21145 


120.5 C 


120.1 


5P3/2 


20874 


20902 






4d 3/2 


14336 


14362 






4d 5/2 


14335 


14360 






6«l/2 


13558 


13509 


239.18(3)"* 


239.2 


6P1/2 


9976 


9973 


39.11(3) e 


38.87 


6P3/2 


9898 


9894 







"Ref. 2C|. 
''Ref. 21 j. 
c Ref. 2j|. 
d Ref. 2S | . 
e Ref. [24 1 . 



TABLE II: Electric dipole transition amplitudes (reduced ma- 
trix elements in a.u.) for low states of Rb. 



Transitions 


Exper." 


Calc. 


5S1/2 - 5pi/ 2 


4.231(3) 


4.246 


5s 1/2 - 5p 3/2 


5.977(4) 


5.994 



a Ref. \25 



Here N is the number of neutrons, Z is the number of 
protons. 

To calculate the PNC amplitude we use the methods 
developed in our previous works 0, |26| . The all-order 
correlation potential £ 27] is used to construct the so- 
called Brueckner orbitals (BO) for the external electron. 
BO are found by solving the Hartree-Fock-like equations 
with an extra operator S: 



(H 



(BO) 



0. 
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Here Hq is the relativistic Hartree-Fock Hamiltonian, 
index a numerate valence states. The BO V , i B °'' and 
energy e a include dominating higher-order correlations. 
The parity non-conserving weak interaction as well as 
the electric dipole interaction of the atom with laser light 
are included in the framework of the time-dependent 
Hartree-Fock approximation 26] which is equivalent to 
the well-known random-phase approximation (RPA). 

In the RPA method, a single-electron wave function in 
external weak and El fields is 



ij) = Vo + Sip + Xe- luJt + Ye luJt + SXt 



- 5Ye iult 



, (4) 



where ?po is the unperturbed state, Sip is the correction 
due to weak interaction acting alone, X and Y are correc- 
tions due to the photon field acting alone, SX and SY are 
corrections due to both fields acting simultaneously, and 
lo is the frequency of the PNC transition. The correc- 
tions are found by solving the system of RPA equations 



self-consistently for the core states 

(H - e c )Sip c = -(H w + SV w )ipoc, 

(Ho - e c - lu)X c = -(H E i + 6V E i)ip c, 

(Ho -e c + u)Y c = -(H E1 + SV^oc, (5) 

(Hq - e c - uj)5X c = -SVsiSipc - SV W X C - SVeiw^Qc 

(Ho -e c + Lo)6Y c = -6V E1 Sip c - SV W Y C - 8V Elw ipo c , 

where index c numerates core states, Hyy is either Hsi 
or Hsd (sec Eq. []}, SVw and SVei are corrections to 
the core potential due to the weak and El interactions, 
respectively, and SVeiw is the correction to the core po- 
tential due to the simultaneous action of the weak field 
and the electric field of the photon. 

The PNC amplitude between valence states a and b in 
the RPA approximation is given by 

^PNC = (lpb\H E l + SV E l\Slpa) + 

(ij b \H w + 5V w \X a ) + ^ b \5V Elw \i) a ) = 
(A\HEi+SV B1 \Sip a )+ (6) 
(Sip b \H E1 + 5V E i\i> a ) + (i>b\SV Elw \ip a ). 

To include correlations in the calculation of the PNC 
amplitude one needs to use BO for the valence states a 
and b in ([6]). The corrections 5ip a and 5ipb to BO a and 
b are also found with the use of the correlation potential 
t: 



(H 



+ Z)6ip a = -(H w + SV)lPoa- 



(7) 



Note that the correlation potential S is the energy- 
dependent operator. To calculate a BO and corrections 
to it one should use the correlation potential at the en- 
ergy of this state, i.e. S = S(e a ) in (J3j) and (JTJ). 

The way of calculation of the PNC amplitude described 
above does not involve direct calculation of the electric 
dipole transition amplitudes or weak matrix elements or 
even the energies, apart from the energies of the 5s and 
6s states. However, it is instructive to make compar- 
isons with available experimental data to have an idea 
of the accuracy of the calculations. For this purpose 
we have performed the calculations of the energies and 
magnetic dipole hyperfine structure constants of the low- 
est s and pi/2 states of Rb as well as the electric dipole 
transition amplitudes between these states. The calcu- 
lations are done with the use of the same approach and 
the same all-order operator S as for the PNC calcula- 
tions [28| , The results for the energies and the hyperfine 
structure are presented in Table U for El transition am- 
plitudes in Table [TXJ Comparison with experimental data 
shows that the accuracy of the calculations is about 0.1% 
for the energies, 0.4 - 0.6% for the hfs and about 0.3% 
for the El transition amplitudes. If we assume that the 
square root of the product of hfs constants of s and p 
states can be used as a test for the weak matrix elements 
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TABLE III: Contributions to the parity non-conserving elec- 
tric dipole transition 5s — 6s in 87 Rb and 6s — 7s in 133 Cs 
(W- 12 iea B (-Qw/N)). 



V_yLJllLI 1ULI L1UI1 


Rb 

a.u. 


% 


wo 

a.u. 


% 


RPA 


1.345 


97% 


8.899 


99% 


Correlations 


0.054 


4% 


0.173 


2% 


Subtotal 


1.400 


101% 


9.072 


101% 


Breit 


-0.006 


-0.4% 


-0.055 


-0.6% 


QED 


-0.003 


-0.2% 


-0.029 


-0.3% 


Neutron skin 


-0.0008 


-0.06% 


-0.018 


-0.2% 


Total 


1.390 


100% 


8.970 


100% 



((s|W|p) ~ y/A s A p ) then the accuracy for the weak ma- 
trix elements is also on the level of 0.3%. Note however 
that the accuracy of this test is limited. For example, 
the value of the ratio (s\W\p) / ^ A S A P is 4% different 
in Hartree-Fock and RPA approximations. This is be- 
cause core polarization effects are significantly different 
for weak and hfs interactions. Only s and P1/2 states 
contribute to the core polarization for the weak matrix 
elements. In the case of hfs interaction the p 3 / 2 states 
also give a significant contribution. Since weak matrix 
elements are simpler, the accuracy for them is expected 
to be higher than for the hyperfinc structure. 



RESULTS AND DISCUSSION 

The value of the spin-independent PNC amplitude for 
the 5s - 6s transition in 87 Rb (without Breit, QED and 
neutron skin corrections) is 



\E 



PNC I 



(8) 



= 1.400 x 10- LZ ea B (-Qw/N). 
This is in very good agreement with the value 
|£pnc| = 1-39(2) x lO- 12 ea B (-Qw/N). 



presented in our early calculations [10j. 

Below we will discuss and compare different contribu- 
tions to the spin-independent PNC amplitudes in rubid- 
ium and cesium and point to some advantages of using 
rubidium in searching for new physics beyond the stan- 
dard model. The contributions are presented in Table HTT1 
We use the 87 Rb isotope as an example. 

Correlations. The total correlation correction to the 
PNC amplitude is small for both atoms. It is 2% for 
cesium and 4% for rubidium. The small value of the 
correlation correction is the result of strong cancellation 
between different terms. This is illustrated by the data in 
Table IIVI where correlation corrections are presented for 
each term in ((6]) . In this table we use notation d = Hex + 
SVei for short. The largest correlation corrections are for 
those terms which have Sip weak correction to the ground 
state. Corresponding corrections for Rb are larger than 



TABLE IV: Correlation corrections (A) to the PNC am- 
plitude, comparison between rubidium and cesium. Units: 
10 _12 ieas(— Qw/N). Indexes a and b stand for 7s and 6s for 
Cs and for 6s and 5s for Rb. 



Approxi- 


{8ip a \d\ip b ) 


(V>o \d\5ipb) 


(lpa\SVBlw\ll>b) 


Total 


mation 














Cesium 






RPA 


-3.041 


11.965 


-0.0249 


8.899 


BO° 


-3.358 


12.454 


-0.0242 


9.072 


A 


-0.316 


0.489 


-0.001 


0.173 


A (%) 


-3.5% 


5.4% 
Rubidium 


0.0% 


2.0% 


RPA 


-0.408 


1.756 


-0.003 


1.345 


BO a 


-0.463 


1.866 


-0.003 


1.400 


A 


-0.055 


0.011 


0.000 


0.055 


A (%) 


-3.9% 


7.8% 


0.0% 


4.0% 



"Brueckner orbitals with core polarization. 



for Cs. Indeed, the closer the valence electron to the 
core the larger the correlation correction. The ionization 
potential for Rb is larger than for Cs. This means that 
the valence electron in Rb is closer to the core than in 
Cs. Also, the cancellation between the contributions of 
6ip a and 5ipb in Rb is not as strong as in Cs. 

Note that the strong cancellations between correlation 
corrections to different terms in (|6]) do not mean poor 
numerical accuracy. The theoretical uncertainty of the 
PNC amplitudes is mostly due to missed terms while the 
numerical accuracy for all terms in ([6]) is high. 

Table IIVI does not include non-Brueckner correlation 
corrections, such as structure radiation, weak correlation 
potential and renormalization of the wave functions [6]. 
These corrections are suppressed by a small parameter 

^valence IE CO 1 



1/10 where E valence a,i\& E core 



are typi- 



cal excitation energies of the valence and core electrons. 
Moreover, their total contribution is practically zero for 
Cs due to cancellation between different terms. It is 
expected to be very small for Rb as well since all other 
relative contributions to Ep^c in Rb and Cs are simi- 
lar. These terms can be calculated if progress is made 
with the measurements. At the moment we just assume 
that they do not contribute to the PNC amplitude or its 
uncertainty. 

The data in Table HVl show that the correlation correc- 
tion to the PNC amplitudes are similar in cesium and ru- 
bidium. Therefore, similar uncertainty is expected. The 
uncertainty for cesium is 0.4 - 0.5% p, 0|- It is natural 
to expect the same uncertainty for rubidium. 

Breit interaction. In was demonstrated in (29j that 
Breit interaction gives significant contribution to the 
PNC in many-electron atoms. The contribution of Breit 
interaction is about -0.6% to the PNC amplitude in 
Cs [H E3] (see also Table HIIj>. To calculate Breit cor- 
rection in Rb we use the same approach as in our pre- 
vious works [30, 311. The Breit Hamiltonian includes 
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magnetic and retardation terms. The Coulomb interac- 
tion everywhere in the calculations is replaced by the 
sum of Coulomb and Breit terms V — > Vq + Vb • Second- 
order correlation correction operator £ is used to cal- 
culate Brueckner orbitals. The correction is found as a 
difference between two results for Eq. ©, one with the 
Breit interaction included and another when it is not in- 
cluded. The resulting Breit correction is about -0.4% of 
the PNC amplitude in Rb. Its relative value is about 1.5 
times smaller than in Cs. Therefore, the uncertainty as- 
sociated with this correction is also smaller for Rb than 
for Cs. 

QED corrections. The quantum electrodynamics cor- 
rections to the Epnc for rubidium are calculated using 
the sum-over-states method. The sum which needs to be 
evaluated is 



PNC 



E 



(6s\d\np 1/2 )(np 1/2 \H w \5s) 
(Qs\H w \np 1/2 )(np 1/2 \d\5s) 



Ea 



En 



(9) 



Here tilde means that core polarization is taken into ac- 
count (e.g. d = Hei +SVei)- Correlations are taken into 
account by using Brueckner orbitals for all ns and np\i 2 
states. 

The QED corrections to the PNC amplitude can be 
seen to arise from three different sources. There are the 
corrections to the El dipole matrix elements (d), weak- 
interaction matrix elements (Hyy), and the energy de- 
nominators. Corrections to the weak matrix elements 
have been considered previously 32, 33| (see also [34l|). 



From these works, we determine the QED contribution 
to the PNC amplitude coming from corrections to the 
weak matrix elements elements to be —0.30(2)%. 

For the corrections coming from the energy denomina- 
tors and dipole amplitudes,we use the "radiative poten- 
tial" method proposed in [8|. By calculating the dom- 
inating terms in equation ([9]) both with and without 
QED corrections, we determine the correction coming 
from the energy denominators to be —0.25% and from 
the dipole amplitudes to be +0.31% giving a combined 
shift of +0.06(3)% for dipoles and energies. 

Therefore we find the total QED shift to the 5s-6s PNC 
amplitude in rubidium to be —0.24(4)%. As expected, 
the Rb result is smaller than that in Cs atom. More im- 
portantly, omitted higher order corrections in Za should 
be much smaller in Rb. 

Neutron skin. The neutron skin correction to the 
PNC amplitude is due to the fact that nuclear density 
in the weak interaction Hamiltonian (fTJ) is not the same 
as nuclear charge distribution. This density is dominated 
by neutrons and if the neutron distribution radius differs 
from the radius of the proton distribution this would lead 
to a correction to the PNC amplitude. It was found from 



TABLE V: PNC 

\5s,Fi) -> \6s,F 2 ) 
10 _11 jea . 



amplitudes (z-components) for the 
transitions in 85 Rb and 87 Rb. Units: 



Isotope 




I 


Fi 


F 2 


PNC amplitude 


Bb Rb 


-44.76 


2.5 


2 


2 


0.091(1 + 0.153*) 








2 


3 


-0.102(1 - 0.206*) 








3 


2 


-0.102(1 + 0.250*) 








3 


3 


-0.137(1 - 0.109*) 


87 Rb 


-46.74 


1.5 


1 


1 


0.071(1 +0.105*) 








1 


2 


-0.123(1 - 0.125*) 








2 


1 


-0.123(1 + 0.167*) 








2 


2 


-0.143(1 - 0.063*) 



an analysis of data for antiprotonic atoms [35| that the 
root mean square radius of the proton and neutron dis- 
tributions differ by 



Ar„ p = (-0.04 ± 0.03) + (1.01 ± 0.15) 



N - 



A 



fm. (10) 



Using this data to correct nuclear density in (TTJ) and re- 
calculating the PNC amplitude leads to -0.2% correction 
for Cs and -0.06% correction for Rb (see Table [TTT]) . Here 
again the correction is much smaller for Rb than for Cs 
leading to smaller uncertainty in the PNC amplitude. 

Combining all corrections we obtain the final value of 
the 5s - 6s nuclear-spin-independent PNC amplitude in 
87 Rb: 

|£+ NC | = 1.390(7) x l0- 12 ea B (-Q w /N). (11) 

and for 85 Rb: 

|£pnc| = 1.333(7) x 10- 12 ea B (-Q w /N). (12) 

We assume the 0.5% uncertainty as it has been discussed 
above. 

The nuclear spin-dependent PNC amplitudes for tran- 
sitions between different hyperfine structure components 
of the 5s and 6s states of 85 Rb and 87 Rb are presented 
in Table [V] 
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